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ON GAUSSIAN BEAMS DESCRIBED BY JACOBI'S EQUATION 

STEVEN T. SMITH 

Abstract. Gaussian beams describe the amplitude and phase of rays and are widely used 
to model acoustic propagation. This paper describes four new results in the theory of Gaussian 
beams. (1) It is shown that the Gerveny equations for the amplitude and phase are equivalent to the 
classical Jacobi Equation of differential geometry. The Gerveny equations describe Gaussian beams 
using Hamilton-Jacobi theory, whereas the Jacobi Equation expresses how Gaussian and Riemannian 
curvature determine geodesic flow on a Riemannian manifold. Thus the paper makes a fundamental 
connection between Gaussian beams and an acoustic channel's so-called intrinsic Gaussian curvature 
from differential geometry. (2) A new formula tt{c/c")^'^ for the distance between convergence 
zones is derived and applied to the Munk and other well-known profiles. (3) A class of "model 
spaces" are introduced that connect the acoustics of ducting/divergence zones with the channel's 
Gaussian curvature K = cc" — {c')-^. The model SSPs yield constant Gaussian curvature in which 
the geometry of ducts corresponds to great circles on a sphere and convergence zones correspond 
to antipodes. The distance between caustics 7r(c/c")^' ^ is equated with an ideal hyperbolic cosine 
SSP duct. (4) An intrinsic version of Cerveny's formulae for the amplitude and phase of Gaussian 
beams is derived that does not depend on an extrinsic, arbitrary choice of coordinates such as 
range and depth. Direct comparisons are made between the computational frameworks used by the 
three different approaches to Gaussian beams: Snell's law, the extrinsic Frenet-Serret formulae, and 
the intrinsic Jacobi methods presented here. The relationship of Gaussian beams to Riemannian 
curvature is explained with an overview of the modern covariant geometric methods that provide a 
general framework for application to other special cases. 

Key ^vords. Paraxial ray, Gaussian beam, acoustic ray, Jacobi equation, Gaussian curvature, 
Riemannian curvature, Hamilton-Jacobi equation 
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1. Introduction. This paper uses Jacobi's equation to derive new formulae for 
the geometric spreading loss and phase through Gaussian beams, and thus provides 
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an alternate method for paraxial ray tracing. [HI [9l [TBI [13 [El [221 [H] The new for- 
mulation, though mathematically equivalent to well-known expressions, provides new 
geometric insight into the physical and intrinsic geometric characteristics of Gaussian 
beams and their relationship to the Gaussian and Riemannian curvature of the prop- 
agation medium. Thus the expressions introduced in this paper establish the connec- 
tion between Gaussian beams and the intrinsic geometry of the propagation medium. 
Four new geometrically-motivated ideas for ducting are presented: (1) the Cerveny 
equations for the amplitude and phase of Gaussian beams are shown to be equivalent 
to the Jacobi Equation of differential geometry; (2) a new formula 7r(c/c")^/^ for the 
distance between convergence zones is derived, where c(z) is the sound speed profile 
(SSP); (3) the intrinsic geometry of acoustic ducting is shown to be equivalent to 
great circles on a sphere with convergence zones corresponding to antipodes; (4) a 
coordinate-free "intrinsic" version of Cerveny's formulae for the amplitude and phase 
of Gaussian beams is presented. 

Paraxial ray methods are generally known as "Gaussian beams" because each 
ray is treated as representing a volume or ray tube in which the ray's amplitude and 
phase in the transverse plane perpendicular to the ray's tangent is determined by 
a Gaussian density. Longitudinally along the ray, Jacobi's equation determines the 
geometric spreading loss, expressed using Riemannian or sectional curvatures, or, in 
the case of 2-d rays, the Gaussian curvature. Therefore, the name "Gaussian beam" 
is highly suitable because Gaussian beams are completely determined by Gauss's 
eponymous curvature. 

One interesting example of a new physical insight derived from Jacobi's equation 
is a simple formula for the distance between caustics: it is shown that the half- 
wavelength distance is about 7r(c/c")^/^, a quantity that depends entirely on the SSP. 
It will be proved that the distance between convergence zones for a Munk profile 
with parameter e and scaled depth W meters is about TTe~^^'^W meters. For SSPs 
with an idealized hyperbohc cosine profile c(z) = cq cosh(z — zo)/W, this distance is 
shown to equal exactly ttW meters for all rays. Caustics arise with positive Gaussian 
curvature; when the curvature is negative, rays diverge and Jacobi's equation quanti- 
fies their divergence, or transmission loss. For example, for linear SSPs with slope c' 
the geometric spreading at time i^ c'^^ is about ct + ^c{c')'^t^. These results are 
both a direct consequence of the fact that the Gaussian curvature of the propagation 
medium equals K = cc" — (c')^. Classifying the SSP by its Gaussian curvature will 
allow for the introduction of model spaces for convergent ducts whose curvature is 
constant positive, divergence zones with constant negative curvature, and simple non- 
refractive spreading with vanishing curvature. Another feature arising from this work 
is an accounting of the additional spreading loss for either reflected or transmitted 
rays at an interface, at which point the Gaussian curvature is not defined. 

Direct comparisons are made between the computational frameworks derived from 
the three different approaches to Gaussian beams: (1) Snell's law, [H [23 HI HI [HI 
[T5l [30] (2) a variant of the extrinsic Frenet-Serret established by Cerveny and col- 
leagues, [H [9] and (3) the new intrinsic methods presented here. Bergman, [H [21 [3] 
apparently the first to recognize the application of Jacobi's equation to ray tracing, 
recently adopted methods from General Relativity to address the problem of comput- 
ing ray amplitudes in a relativistic acoustic field. The non-relativistic intrinsic results 
developed in this paper are equivalent to many of Bergman's if one uses the space-like 
part of his pseudo- Riemannian Lorentz metric. 

It is perhaps noteworthy that Jacobi's equation and its full implications for Gaus- 
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sian beams has, apparently, not yet appeared in the acoustical literature. This lacuna 
might be defended by a few historical observations. Lord Rayleigh was initially dis- 
missive of the practical applications for acoustic refraction, noting in 1877 "almost the 
only instance of acoustical refraction, which has practical interest, is the deviation of 
sonorous rays from rectilinear course due to the heterogeneity of the atmosphere." [24] 
Though the foundations of non-Euclidean geometry had been established by this time, 
differential and Ricmannian geometry remains relatively little known in engineering 
applications to this day, in spite of the fact that covariant analysis is the natural ap- 
proach to many physical and engineering problems, as will be seen in its application 
to Gaussian beams here. 

In Section [2] the ray equations developed using the classical Euler-Lagrange, 
Frenet-Serret, and Hamilton- Jacobi formulations, followed by the amplitude and phase 
along a specific ray using Cerveny's Hamilton- Jacobi approach. This standard devel- 
opment is then recast using an intrinsic parameterization that will be shown equivalent 
to Jacobi's equation. Section [3] develops an intrinsic formulation of Gaussian beams 
based on Jacobi's equation and explores the physical consequences of this approach, 
including a computation of the distance between convergence zones and a classification 
of SSPs using "model spaces" of constant Gaussian curvature. 

2. Gaussian Beams in Horizontally Stratified Isotropic Media. The sim- 
plest and most frequently encountered case of acoustic rays in a horizontally stratified 
isotropic medium will be analyzed first, both because of its importance to applica- 
tions, as well as to fix ideas and provide comparisons to widely known results. As 
usual, denote the three spatial coordinates by the variables x, y, and z, time by t, the 
spatial infinitesimal arclength by ds^ = dx^ + dtp -\- dz^ , and the depth-dependent SSP 
by c{z). Assume that c(z) is differentiable except at a discrete set of points where ei- 
ther c{z) itself is discontinuous (e.g., Snell's law), or its first derivative c'(z) — dc/dz 
is discontinuous (e.g., method of images at a boundary). These special cases are 
treated separately below. It will be convenient to use the cylindrical coordinates 
r — yjx'^ -I- y^ (range), = tan^^(y/a;) (azimuth), and z (depth, pointing down), so 
that ds^ = dr^ -\- r^ dcfi^ -t- dz^ . Let x denote the 3-by-l vector of chosen coordinates 
(we will use cylindrical coordinates here, but the physical solution in independent of 
this choice, and we will generalize to arbitrary choices below). The time required to 
travel along an arbitrary continuous path x(i) — (r{t), (f>{t), z{t)) equals 

T[x(t)] =Jdt = J-^^^dt = J c~\z)ir^ + r^^ + z')'/' dt. (2.1) 

where x(i) = dx./dt and ds/dt = c{z). 

To treat the general case of geometric spreading loss, the following standard 
notation for metrics will be used. The quadratic function 



whose square-root appears in the functional integral of Eq. (2.1) is a Ricmannian 
metric representing the square infinitesimal intrinsic path length. This so-called Fer- 
mat metric [with coefficients G = {gij)] induces an inner product (x,y) = ^(x, y) = 
c~^(z)(fif2 -f-r'^0102 + iii2) and norm ||x|p = 5(x,x) on the space of tangent vectors 
(r,0, i)'^ at each point (r, 0, z)'^. Note that rays arise from time-minimizing paths 
w.r.t. the Fermat metric, whereas straight lines are extremals for the standard ex- 
trinsic inner product (x, y)2 = Tli^^V^ ^'^'^ 2-norm ||x||2 = 'Ylii^^Y that do not have 
position-dependent speed. 
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Fig. 2.1. Rays for the Munk SSP ilSj 
with with parameter e = 0.00737, scaled depth 
z = {z — zc))/W with W = 650 m, and (noncon- 
stant) Gaussian curvature K = ec^/W^ at zq, 
yielding convergence zones at about every 
TTe~^'^W = 23.8km as illustrated by the red 
dots (»). The rays are computed using a stan- 
dard 4th order Runge-Kutta ODE solver. 
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Fig. 2.2. Rays for the hyperbolic co- 
sine SSP c{z) = Co cosh(2 — 2o)/iy with zq = 
W = 1 km. This SSP yields a space of con- 
stant positive curvature, resulting in a duct 
with convergent rays and caustics independent 
of any initial elevation angle Bq at ranges 
ttW = 3.14159 ...km (Theorem\3^, illustrated 
by the red dots (•). The rays are computed using 
a standard J^th order Runge-Kutta ODE solver. 



Before the main results of the paper involving the geometry of acoustic ducting 
and divergence are presented, the standard ray equations must be derived. Subsec- 
tions |2.1| and |2.2| derive the ray equations using the Euler-Lagrange, Frenet-Serret, 



and Hamilton-Jacobi formulations, then subsection 2.3 derives the amplitude and 
phase along a spe cific ray using Cerveny's Hamilton-Jacobi approach. Concluding 
this section in 



2.4 



Cerveny's equations will be recast using an intrinsic parameteriza- 
tion that will be shown in section |3] to yield Jacobi's equation expressed in form that 
yields geometric insight into acoustic ducting. 

2.1. Ray Equations: Euler-Lagrange and Prenet-Serret. Huygens's prin- 
ciple and Fermat's principle of least time implies that rays satisfy the Euler-Lagrange 
equation 



d dL 
dt 9x 



9x 



0. 



(2.3) 



In the case of freely propagating rays, the initial conditions x(0) = xo and x(0) = 
xg are specified. In the case of eigenrays, the boundary conditions x(0) = xq and 
x(end) = Xcnd are specified. The radial symmetry in the case of a horizontally 



stratified medium along with Eq. (2.3) implies that in cylindrical coordinates, (j){t) = 
00 J so we will restrict our attention to the (r, 2;)-plane, as usual. Applying Eq. ( |2.3[ ) 
to the Lagrangian 



L(t; r, z; f, i) = c-^{z){f'^ + z^fl"^ 



(2.4) 



of Eq. (2.1 ) such that (/) = 0, rays in a horizontally stratified medium satisfy the well- 



known differential ray (Christoffcl) equations 



f — 2(c /c)fz = 0, 
z + {dlc){r'^-P)=0. 



(2.5) 
(2.6) 



See Figs. |2.1| and |2.2| for computed rays determined by the Munk and hyperbolic cosine 
SSPs. Parameterization of rays by travel-time t is said to be "natural" or "intrinsic" 
because rays minimize travel-time. However, it is oftentimes computationally conve- 
nient to parameterize rays by "extrinsic" arclength ds = {dr"^ 4- dz'^)^''^ , in which case 
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Eqs. (2.5)-(2.6) become the first Frenet-Serret formula 

(P /r\ _ c^dr / dz/ds \ 
ds^\z)~'cds \-dr/ds) ' 

i.e., dT/ds = >cN, wlicre T = {dr/ds, dz/ds) is tire ray's tangent vector, 
(dz/ds, —dr/ds) is its normal (always defined in the same direction), and 

K = [d / c){dr / ds) = — Cn/c 

is the ray's extrinsic curvature, and 

c„ = (9c/ax)-N = -cJ (dr/ds) 



(2.7) 



N 



(2.; 



(2.9) 



is the first derivative of the SSP c in the normal direction. Physically, rays travel 
along time-minimizing "straight" lines, and our arbitrary choice of describing th e ray 
using the coordinates (r{s),z{s)) and the resulting extrinsic curvature of Eq. (2.81 



has no physical meaning to the ray itself; different coordinates would yield a different 
curvature. This traditional extrinsic approach will be compared below with an intrin- 
sic approach that uses the intrinsic Gaussian curvature of the propagation medium. 
The quadratic coefficients for the tangent vector (first derivative) terms that appear 
in Eqs. (2.5)-(2.7) are called Christoffel symbols (of the second kind), and are crucial 
in quantifying the amplitude and phase along the ray caused by geometric spreading. 
U sing t he H amiltonian formulation developed in the next subsection, the ray equa- 
tions ([2!^-([2!6| for the initial conditions (r(0),z(0)) = (0,zo)^ and (f(0),i(0)) = 
c(zo)(cos6'o, — sin^o)^ can be reduced to a function r{z;9o) expressed using a single 
integral, which represents a continuous version of Snell's law: [H [6l [TTl [151 130] 

ac{z,) 



r{z;eo) 



1 



^1 1^1 



(^.)) 



1/2 



dzi, 



(2.10) 



where a — c ^ 



Zo)cos0o is the Snell invariant. It is a useful exercise to show that 
Snell's law is a consequence of Noether's theorem. 

2.2. Ray Equations: Hamilton-Jacobi Equation. The Euler-Lagrange equa- 



tions (2.5)-(2.6) or (2.7) determine the ray itself, and the Hamiltonian or canonical 
form of these equations determines the amplitude and phase along the ray. Pa- 
rameterizing the ray by the range r, corresponding to the Lagrangian L(r; z; z') = 
c-i(z)(l-F0'2)i/2 with z' = dz/dr, the time of travel T[r; z] = /J' c-^{z){l + z'^)^/'^ dr, 
along a ray from (0, zq)'^ to the point (r, z)'^ satisfies the Hamilton-Jacobi nonlinear 
partial differential equation 

dr V ' ' 9z 



0, 



(2.11) 



where H{r; z; Q = (z' — L{r; z; z') and C = dL/dz'. Rays satisfy the canonical 
equations z' = dH/dC,, C,' = —dH/dz, and Hamilton-Jacobi theory provides the rela- 
tionship c 



dT/dz along rays. For the Lagrangian of Eq. (2.4 1 

^2\l/2 



Hir;z:C) = -{c-'iz)-CT^ (2-12) 

and the Hamilton-Jacobi equation of Eq. ( |2.11 1 is simply the well-known eikonal 
equation 



dT 

dr 



^ (dT 



\^\ 



(2.13) 
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Fig. 2.3. Ray tube in a horizontally stratified medium. The azimuthal component, not shown, 
is perpendicular to the page. All distances in the figure represent actual, extrinsic Euclidean dis- 
tance. The extrinsic distance between nearby rays along the normal vector N is denoted as 5q, 
and the nearby ray's additional length equals ■;^c{5p/5q) &q^ = i^c(p/q) Sq^ , where q = WvciSq/SBo is 
the extrinsic geometric spreading, and 5p = dT/dSq and p = dT/dq are the conjugate momenta 
corresponding to &q and q. 



2.3. Paraxial Ray Equations: Hamilton-Jacobi Form. The amplitude along 
rays is determined by the spreading of nearby rays, and the phase away from the ray 
is determined by the time difference. Thus, rather than the skeletal objects deter- 



mined by Eqs. (2.5)-(2.6l or (2.71, rays are viewed as tubes or beams possessing both 
amplitude and phase. As first established by Cerveny and colleagues, [51 1^ [TCI E^ [T5] 
the equations for the amplitude and phase along a ray tube are given by the canonical 
equations after a convenient change of variables involving the ray itself. Let s be the 
arclength along the ray and let 5q be the small or infinitesimal distance away from 
the ray, measured perpendicularly along the normal N at arclength s such that 



r{s,5q) 
z{s,6q) 



r{s) 
z{s) 



Sq 



dz/ds 
~dr/ds 



(2.14) 



and let 



5p = dT/d6q 



(2.15) 



be the small or infinitesimal derivative of travel-time w.r.t. 6q from the ray along the 



line N, all illustrated in Fig. 2.3 Thus, 5q quantifies the spread of nearby rays, and 
therefore determines the ray's amplitude, and as shown by Cerveny [^ the conjugate 
momentum 6p appears in the difference in time-of-travel nearby rays along N and thus 
determines the ray's phase. Following Cerveny [5] (also see Wolf and Krotzsch |31) ) 
Hamilton-Jacobi theory is used to determine the governing equations for Sq and 5p. 
By the chain rule dT/d{s,Sq) = dT/d{r,z)-J, where J = d{r,z)/d{s,Sq) = 
(/iT, N) is the Jacobian matrix along the ray at distance s at Sq — and h is the 
scale factor 



h{s, Sq) = 1- [{d / c){dr / ds)) \sq=o Sq = l + {cn/c)\&q=o Sq, 



(2.16) 



the left-hand-side of the eikonal equation Eq. (2.13) may be computed in the ray- 
centered coordinates s and Sq: 



1 fdT 
h^{s,Sq) \1^ 



dT 
dSq 



1 

c'^{s,Sq)' 



(2.17) 
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By Hamilton-Jacobi theory, the Hamiltonian describing these paths is, 

H{s; Sq; 5p) = -/i(s, Sq){c-\s, Sq) - Sp'^^^^ 
For small 6q and Sp up to second-order, 



c{s,dq) Cs cl 2 V cl cf 



H{s- Sq; Sp) ^ ^hc~\s, 5q){l - \cl Sp^) + ■■■ 



2c2 



5q' 



:CsSp 



where 



c(z(s)j, and 



c"{z{s)){dr/ds 



(2.18) 

(2.19) 
(2.20) 
(2.21) 

(2.22) 

-dH/dSq 

to Eq. ( 2.21 1 yields the paraxial ray tracing equations, a system of ordinary differential 
equatioiip El [H [16l [H [HI EH [25] 

^=c,Sp, SqiO) = 0, (2.23) 

as 

dSp 

~ds ~~ 2cl 

the initial conditions arising from the facts that the spread of rays separated by a 
small angle S9[) is zero at s = 0, and the differential time equals ST — Cq SqSO^. 

2.4. Paraxial Ray Equations: Sturm-Liouville Form. Before explaining 
the full significance of the complete second-order terms in the paraxial ray equations 



is the second derivative of the SSP in the direction of the ray normal. 
Applying Hamilton's equations of motion dSq/ds = dH/dSp, dSp/ds 



Sq, 



Sp{O)^CoU0o, 



(2.24) 



in the next section, it will be helpful to express Eqs. ( 2.23 )-( 2.24) as a single second 



order differential equation using the natural or intrinsic parameters of the problem. 
Consistent with much of existing literature on paraxial rays, rays are parameterize 
using the arclength s. However, physically rays do not minimize arclength, but travel- 
time i, and therefore the paraxial ray's underlying physical and geometric properties 



will be revealed by using the intrinsic parameterization dt 



ds. Also the variables 



Sq and Sp = dT/dSq are expressed using arclength or distance; the corresponding 
intrinsic variables arise by converting distance Sq to travel-time Sq and slowness Sp to 
the dimensionless Sp using the sound speed c: 



Sq = c Sq, and 



f ~ dcf f. 

Sp = cSp. 



(2.25) 



Note that Sq is the travel-time between nearby rays and Sp — dT/dSq is the travel- 
time derivative along a straight-line path normal. The geometric significance of these 
intrinsic coordinates is made clear by the following new theorem: 

Theorem 2.1. Let Sq and Sp satisfy Eqs. {2.23)-(2.24-). Then these extrinsic 
first-order paraxial ray equations are equivalent to the intrinsic second-order paraxial 
ray equation 



Sq+{cc" -{c'f)Sq^Q 



(2.26) 



where Sq ~ c ^ Sq, Sp = cSp, Sq = (d^ /dt^)Sq, and the coefficient cc" — (c')^ is 
the Gaussian curvature of the propagation medium w.r.t. to the Fermat metric of 
Eq. K§. 
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The proof involves parameterizing by travel-time t and expressing the Hamilton 



form of the paraxial ray equations of Eqs. ( 2.23 )-( 2.24) as the coupled first-order 

system, 

6q = -c'{dz/ds) 6q + 6p dq{0) = 0, (2.27) 

dp = -ccnn Sq + c!{dz/ds) 5p 5p{Q) = SOq. (2.28) 

Taking another derivative of the first equation yields 6q + K 6q = with K = cc" — (c')^- 

The theorem is thus proven by demonstrating that this expression for K is pre- 



cisely the Gaussian curvature w.r.t. the Fermat metric of Eq. (2.2). Given an ar- 
bitrary metric 5 on a 2-dimensional manifold with coefficients gij = g{X.i, Xj) w.r.t. 
a basis { Xi }, the Gaussian curvarure is determined |26j by the expression K = 
-Ri2i2/(5ii522 - gh) in which i?„H = Y.n9inR"'jki, R'^jki = {drij/dx'')-{dTlj/dx^) + 
^ni^u^kn ~ ^kj^ln) ^'"^ ^hc coefRcients of the Riemannian curvature tensor R, and 
the Christoffel symbols (of the second kind) T^j — ^J2i 9^^ {{dgu/dx^) + (dg^i/dx'^) — 
{dgij/dx'')) appear in the quadratic coefficients in the ray (geodesic) equations 

i-'^ + E^jT^i^i^ =0 (2-29) 

with (f?*^') = {gij)~^ representing the matrix inverse of the matrix of metric coeffi- 
cients. [11 Hlj The Christoffel symbols 



c'/c or ^ '^' I -c'/c 



can be read directly from the ray equations of Eqs. (2.5|-(2.6|; these along with 
the matrix {gij) = c~^( g -^ 1 can be used directly with the preceding equations to 



compute the Gaussian curvature K — cc" — (c') . As will be proved, Eq. (2.261 is 



precisely the Jacobi equation; therefore, the Hamiltonian form of the paraxial ray 



equations ( 2.23 )-( 2.24) is mathematically and physically equivalent to the intrinsic 
Jacobi equation. This equivalence between Cerveny's formulation of Gaussian beams 
and Jacobi's equation is a new result. The functional form of the Gaussian curvature 
will be used to derive the distance between convergence zones and "model spaces" of 
constant curvature for the SSP. 

3. Intrinsic Geometry of Gaussian Beams. 

3.1. Transverse Amplitude of a Gaussian Beam. This section contains an 
analysis of Gaussian beams that use Jacobi's equation to quantify their amplitude 
and phase. All the results obtained so far have been for rays and infinitesimal vari- 
ations around them. In contrast, Gaussian beams have real, finite physical width 
that must be accounted for in the plane transverse to the ray. The physical width 
has different implications for the Gaussian beam's amplitude and phase. The am- 
plitude is determined by both the geometric spreading loss along the ray and the 
initial amplitude distribution in the transverse plane. The phase is determined by 
the differential change in travel-time along the transverse plane. Mathematically, 
the scaled normal vector 5q N represents the infinitesimal change in ray position as 
a function of the take-off angle for a fixed path length s. Gauss's lemma implies 
that this travel-time is fixed to first-order, but says nothing about the second-order 
terms, i.e., x(s; 6*0) -I- 5qN = x(s; ^o + ^^^o) + 0{59l) for ^^o > 0. This ray-front cur- 
vature effect is seen in Euclidean space by holding a ruler next to a circle or cylin- 
der and observing the size of gap — which for a circle of radius R has a gap size of 
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Rsec66 — R — {R/2)59^ + 0{S9'^). The second-order travel-time difference is often 
physically significant because it is of the order of several wavelengths at a broad range 
of frequencies /o and therefore does contribute to the Gaussian beam's phase 27r/oi. 
Because the travel-time difference in the transverse plane is quadratic, the phase of 
Gaussian beam necessarily has a Gaussian distribution. Furthermore, these second- 
order affects are already accounted for in the Riemannian curvature terms determining 
the Gaussian beam's amplitude, so the Gaussian beam's amplitude in the transverse 
plane is determined by the effect of geometric spreading on initial conditions. In this 
subsection we will quantify the impact of the initial conditions on the amplitude and 
the second-order travel-time differences on the phase. Both quantities are determined 
by the ray distance Sq, which will now be formalized for both extrinsic and intrinsic 
geometric spreading. 

Definition 3.1. (Geometric spreading) Leix(s;0o) = {f{s;do),z{s;6f))y be 
a ray parameterized by arclength s with initial conditions x(0; 6*0) = (0, zo)^ o-^d 
{d/ds)x{O]0o) = (cos^Oj ~ sin^o)^ *^ o. horizontally stratified media with SSP c{z). 
The extrinsic geometric spreading q{s) along the ray caused by a change in elevation 
angle 6q is defined to be 

<z(s) = |l(5x(s;0o)/a0o).ll2, (3.1) 

where ||-||2 denotes the standard 2-norm and the standard notation {d/d9o)s means 
partial differentiation w.r.t. Oq while holding arclength s constant. 

Theorem 3.2. (Geometric spreading for horizontally stratified media [TS]) Let 
x(s) = (r(s; 0o)i •2^(5; ^0)) be a ray in a horizontally stratified medium with a twice- 
differentiable SSP c{z) parameterized by arclength s and with initial elevation angle 6q. 
The extrinsic geometric spreading q and the corresponding canonical variable p = 
dT/dq along the ray are determined by 

q — lim 6q/d0o, p — lim Sp/S0o, (3-2) 

and satisfy the Hamilton equations 

dq/ds = cp, q{0) = 0, (3.3) 

dp/ds = -(c„„/c2)g, p{0) = c^\ (3.4) 

The proof is a trivial application of the definition of the small or infinitesimal 



difference 5q defined in Eq. (2.14 1. The theorem follows immediately by dividing Eqs. 



( 2.23 )-( 2.24) by 69q and taking the limit. Thus we have a complete description of 



the extrinsic geometric spreading for Gaussian beams, including the complete second- 



order terms introduced in Eq. (2.24 1. To appreciate the geometric significance of the 
spread, we will define the intrinsic geometric spreading and equate this concept with 
the Jacobi equation first encountered in Theorem |2.1[ expressed in this new theorem: 
Definition 3.3. (Intrinsic geometric spreading for horizontally stratified me- 
dia) Let (r{t;9o), z{t;9o)^ be a ray in a horizontally stratified medium with a twice- 



differentiable SSP c{z) as in Theorem 3.2, but parameterized by path time t. The 



intrinsic geometric spreading f (t) along the ray is defined to be 

^ (50n->O 09q 

where Sq is the travel-time between nearby rays separated by 69q at their initial point 



defined in Eq. {2.25), and q{t) is the extrinsic geometric spreading defined in The- 



orem 3.2 [Intrinsic geometric spreading is denoted by the Arabic letter 'ayn, ('?•', 
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pronounced like the end of 'nine' spoken with a strong AustraUan accent) in recog- 
nition of the mathematician Ibn Sahl who discovered Sneh's law of refraction around 
984. Ibn Sahl used the symbol f to denote the center of a lens [2?].] 

The intrinsic geometric spreading is a consequence of an obvious corollary to 



Theorem 3.2 arising from Eqs. ( 2.27 1-( 2.28) for the intrinsic variables 



q 



c-\ 



lim 6q/69i) = f- , 

56)0^0 ^ 



and 



p = cp = lim Sp/S9o = f - 

SOn — ^0 



c'sine*^. (3.6) 



Note that the path time derivative p — dT/dq w.r.t. distances along the extrinsic path 
normal N is not exactly equal to the derivative f = (d/dt)^ of the intrinsic geometric 
spreading; the geometric reason for this difference will be explained in Section |3.7[ 
We are now ready to introduce Jacobi's equation and thereby compute the ray's 
amplitude determined by the geometric spreading, as well a Gaussian beam's phase 
in the transverse plane perpendicular to the ray's tangent. 

Theorem 3.4. (Intrinsic geometric spreading for horizontally stratified me- 
dia) Let (r{t; 9q), z{t; Oq)) be a ray in a horizontally stratified medium with a twice- 
dijferentiable SSP c{z) param,eterized by path time t and with initial elevation angle 6q. 
The intrinsic geometric spreading Fjt) along the ray satisfies the Sturm- Liouville equa- 
tion 



l+K{t)l^Q- ^(0) = 0, t(0) = l, 



(3.7) 



where 



K 



(cf 



(3.; 



is the acoustic Gaussian curvature of the Fermat metric g(x, x) = c ^ 
Eq. pjp. 



(z){P + P) 



of 



Eq. (3.7) is simply the intrinsic Jacobi equation in two dimensions encountered 



above in Eq. (2.26). The theorem follows immediately by dividing Eq. (2.26) by 56q 



and taking the limit. In general, the Jacobi equation for an arbitrary Riemannian 
manifold is. 



VlV 



R{V,T)T^O, 



(3.9) 



where W^V is the second covariant derivative along the ray's tangent vector T - 
(f, i)'^ of the variation vector V = 6x/S9o — FY along the unit vector Y = —{z, —r-y 
and the Riemannian curvature tensor equals 



R{X, Y)Z = VxVyZ - VyVxZ - V^x.y] 



(3.10) 



where [X, Y] = XY - YX is the Lie bracket. [lUlIIlISS] Note that the orthonormal 
frame (T, Y) is parallel along the ray, i.e., VtT = [Eqs. ([2|5|-([2!6|] and VtY = 
so that VtV = £y -I- ^tY = ^Y . Furthermore, [T, V] = because the vec- 
tor fields T and V are defined as independent partial derivatives of the ray x(i;6'o); 
therefore, VtV = VvT by the property Vx^ — VyX = [X,l^] of covariant dif- 
ferentiaition. Taking another covariant derivative w.r.t. T and computing an inner 
product with Y along with the definition for Gaussian curvature [THl [131 US] 



K = g{R{Y, T)T, Y)l{g{T, T)g{Y, Y) - g{T, Y) 



(3.11) 



yields Jacobi's equation of Eq. (3.7) 
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Jacobi's equation provides a physically appealing description of Gaussian beams; 
indeed, interpretations of Jacobi's equation and Gaussian curvature provide the fol- 
lowing new insights into the problem of ray acoustics. At depths where the SSP c{z) 
is concave, i.e., a sound duct below which rays refract upwards and above which they 
refract down, cc" > 0, c' sa 0, and the acoustic Gaussian curvature is positive, yield- 
ing a sinusoidal behavior with wavelength — the distance between caustics — of about 
TTcK~^ for the geometric spreading, as is expected and encountered with caustics in 



sound duct (Figs. 2.1 2.2 1. Note well that caustics, defined to be points where the 
geometric spreading vanishes (these are called "conjugate points" in the mathematical 
literature), will occur at integer multiples of the range 

half- wavelength ^^1/9 , , 77x1/9 ,„^^^ 

^ w ncK-^/^ « TT c/c" / . (3.12 

distance ^ ' ^ ^ ' 

(the last approximation valid when the SSP's second derivative dominates). At depths 
where the SSP c(z) is convex, i.e., divergent zones below which rays refract downwards 
and above which they refract up, or at depths with a linear SSP, cc" < 0, (c')^ > 0, 
and the acoustic Gaussian curvature is negative, yielding an exponentially growing so- 
lution to geometric spreading with length constant of about {~K)~^^^, and therefore 
a large spreading loss. Constant SSPs imply that the acoustic Gaussian curvature is 
zero, yielding geometric spreading that simply grows linearly with the path length, 
as expected. At interfaces where the SSP or its first derivative are discontinuous, 
such as the bottom or surface where the method of images is used to model reflec- 



tions, Eq. (3.7) can be integrated using the standard modifications, provided below, 
necessitated by the Dirac delta function. These observations will all be formalized in 
Section [3^ below. It also shown below that the known extrinsic formulae for the geo- 
metric spreading loss based on Snell's law and variants of the extrinsic Frenet-Serret 
formulae 8, 9,[22l[T5] satisfy the Jacobi equation. 



Eq. (3.71 may be proven directly via the second variation, but this direct ap- 
proach, which involves a long, messy computation that provides almost no geometric 
insight, motivates introduction of the cleaner and simpler modern covariant approach. 
Nevertheless, direct computation establishes one important result for readers without 
backgrounds in covariant differentiation and Riemannian geometry, so we will provide 
a sketch in this paragraph. Gauss's lemma assures us that nearby rays are perpendic- 
ular to a given ray's tangent vector (r^z)'^] therefore, nearby rays all take the form 
Sx — (^^{t)z69o, —^{t)r 69o) for the geometric spreading function ^ (t) a long the ray. 



It can be shown using the Hessian matrix of the Fermat metric of Eq. (2.2 ) w.r.t. (x, x) 



the travel-time of Eq. (3.141 equals 



and application of the Christoffel Eqs. (2.5|-(2.6), that the second variation [12] of 



S^L[6x] ^ L[x + (5x] - L[x] - SL[Sx] - o{\\5x\\ 



[ °"' ({^{-c-^cc" - (c'f) + (c7c)z)^2 ^ ^2 ^ 2{c'/c)z^[) dt. (3.13) 



Because these second order terms necessarily satisfy the Euler-Lagrange equation for 
nearby extremals x + (5x, the function (5x must satisfy the Euler-Lagrange equation 



for the second variation, known as Jacobi's equation. Applying Eq. (2.3) for the 



unknown function y{t) to the function appearing in the second variation of Eq. (3.13) 



yields the Jacobi equation of Eq. (3.7). This establishes Theorem 3.4 but it does not 



provide any immediate connection to the problem's geometry or intrinsic curvature. 
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nor does it suggest a method to generalize the result to three (or more) dimensions. 
These connections will now be established within the broadly general framework of 
covariant differentiation. 

3.2. Properties of Geometric Spreading. In this section we will review, in 
the context of the results presented above, some simple properties about geometric 
spreading. First, the geometric spreading q{t) quantifies the distance of nearby rays 
from the ray x(i), parameterized by its initial elevation angle Oq. Second, in the trivial 
case of a constant SSP, geometric spreading simply equals range, q(t) = ct (constant 
SSP). Third, and crucially for understanding the results presented in this paper, 
the direction of geometric spreading is perpendicular to the ray's (extrinsic) tangent 
vector, T — dx/ds, and is therefore also perpendicular to the ray's intrinsic tangent 
vector X = {dx/ds){ds/dt) = c{z)T. This orthogonality property, called Gauss's 
lemma because it is among of the first key observations by Gauss in his establishment 
of the intrinsic theory of surfaces, is also a straightforward consequ ence of Euler's 



first variation of the functional L[x(t)] = J^"""^ F{t,x,x)dt [as in Eq. (2.1|] along an 
extremal x(i) with only the initial endpoint fixed. The first variation formula for an 
arbitrary functional variation (5x(i) such that (5x(0) = is 



SL[Sx.] =L[x + 5x] -L[x] 

p \ fiJ^ tend 

= (9F/ax)Miend)-0, (3.14) 
where the integral vanishes because x is assumed to be an extremal and satisfies 



tend /Qp ^ ^ Qp\ dF *" 

\ 9x dt dx J dx 



the Euler-Lagrange equation (2.3 1, and (5x(0) = because the starting point is fixed. 
The first variation must vanish necessarily for arbitrary functions 6x{t); therefore 
(dF/dx) 5x{t end) = 0. In our specific case of minimum travel-time rays, F{t,x,x) = 
c~^(z)(r^ + py/^ and Sx = {dx/d6o)tS6o and Gauss's lemma follows immediately 



from Eq. p.l4|) for the case of acoustic rays: 

= 0, (3.15) 




Note that this basic argument holds for all extremal paths, whether or not they 
are minimizing or pass through a conjugate point a.k.a. caustic. Consequently from 
Gauss's lemma, a fourth property is that the geometric spreading depends on both r- 
and i-components of the ray's tangent vector. Fifth, there is a geometric spreading 
term in the range-azimuthal (r, (/))-plane, which, by symmetry, is trivially dq{t)/d(f>o = 
r{t). This intuitively obvious fact will be proven as an illustration of the covariant 
formalism introduced below. Sixth, and finally, geometric spreading as defined here 
is physically extrinsic (relative to the metric of travel-times between points) because 
pressure transducers integrate power over physical area, which is why the Euclidean or 



2- norm appears in Definition 3.1 Nonetheless, a simple change of scale [multiplication 
by l/c(z)], converts the problem of geometric spreading to an intrinsic framework that 
depends only upon the travel-time between points, i.e., only on rays. 

The transmission loss along a ray tube caused by geometric spreading is given 
by the infinitesimal area at a standard reference range (r^of = 1 m) divided by the 
infinitesimal area of the ray tube at an arbitrary distance t, both scaled by the inverse 
sound speed along the ray [15] : 

TL (geom.) = lim — = (relm^). (3.16) 

S(j>,se^o corq{t) 69 (50 corq 
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6e„ 



y2C{p/q)bq^ 
bq = qbQ^ 



x{t- Q,+bQ^ 

= qbQ^ = ctbQ, 



mAt)) 




Fig. 3.1. Phase of the Gaussian beam along 
the ray's extrinsic normal N = dT/ds. The dif- 
ferential time Stc a distance 5r] along the normal 
is given by the quadratic Stc = 5 ip/l) ^V^ ~ 
2C~^(^/j_ + c' sinS) 5?7^. All distances in the 
figure represent extrinsic Euclidean distances. 



Fig. 3.2. Geometric spreading expressed 
using ray-centered and extrinsic range and depth 
coordinates. All distances in the figure represent 
extrinsic Euclidean distances. 



When the density also varies with depth, Newton's second law must be included: 

PsCs cos 00 
PoCQrq 



TL (geom. 



(re 1 m^ 



(3.17) 



3.3. Transverse Phase of a Gaussian Beam. The differential time St^ per- 
pendicular to the ray (extrinsic) determines the differential phase e^ ^"^ of the Gaussian 
beam, where dip = 2'KfcStc, fc is the center or carrier frequency of the sound prop- 
agating along the ray, and j = %/— T- In terms of the extrinsic geometric spreading 
variables q and p — dT/dq, the differential time a distance Srj along the ray normal 
at {^r{s),z{s)) is given by the formulae p[|T5] 



2q ' 2 




(3.18) 



where Srj = cSp,. This fact is illustrated in Fig. |3.1[ To compare this extrinsic 
formula with the intrinsic results developed in Section |3.7[ it is worthwhile to give a 



brief geometric proof of Eq. (3.18). By inspection, or formally by Gauss's lemma of 



Eq. (3.15) the differential path time along the ray normal is a quadratic function of 



distance up to second order, i.e. 



SU 



\Qc5rf 



(3.19) 



for some quadratic coefficient Qc to be determined. As seen in Fig. 3.1 the slope of 
this quadratic function at 5rj equals 56, the differential change in angle: 

59 = cd{5tc)/d{6rj)\sy^=sq=qS9a = cQcSq = cqQcS9Q. (3.20) 

However, 69 may be computed directly via an inner product along with Eq. (3.3), 

69= {N,{d/ds)x{s;9o + 69o)) = {N,{d/ds){x{s;9o) + q69oN)) =cp69o (3.21) 



Equating 69 from Eqs. ( |3.20| a nd (|3.21| yields Qc ^ p/g, from w hich E q. ( |3.19[ ) be- 
comes the first equality of Eq. (|3.18|). Applying Eq. (|2.27|) to Eq. (|3.18|) yields 



„2P 






d sin 9 



(3.22) 
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and the second equality. 

Summarizing the results of the previous sections, the transverse amplitude and 
phase of a Gaussian beam emanating a distance Srj normal to the ray and with initial 
source level SL and azimuth end elevation angle fan widths S(j)o and SOq is given by 



P{s■,0o;S^i)^ 



'SL-c, 


cos 6*0 


2 

exp 
J27rf, 


corq 

1 

SL-a 

exp 

. ^^ J 



it+- I - + c' sin 



5fi^ 



(3.23) 
(3.24) 



Note that Gaussian beam is an appropriate name for rays with this functional form. 
A Gaussian distribution of amplitude is typically introduced through the initial con- 
ditions for the Gaussian beam. 

3.4. Munk, Linear, and SSPs v^rith Constant Curvature. 

3.4.1. Munk Profile. The distance between convergence zones/cycle distances 
for the Munk profile c{z) = co(l-l-e(z+e~^ — 1)) with z — {z — zq)/W can be computed 
directly via Eq. (3.12). At z = zq, the Gaussian curvature K = cc" — {c'Y = ^Cq/W'^, 
proving the new theorem: 

Theorem 3.5. The cycle distance of the Munk profile approximately equals 



Munk profile 
CZ distance 



-'/'w. 



(3.25) 



E.g., for the nominal parameters e = 0.00737 and W — 650m, this yields a dis- 



tance of 23.8 km (Fig. [2TJ). 

3.4.2. Linear SSPs. The ideas developed to this point are nicely illustrated 
using the example of Gaussian beams in medium with linear SSP c(z) = cq + "/z, in 
which case rays are simply circles in (r, z)-space: J15j 



r{6) = zq^ sec Oq (sin ^q + sin(0 — ^q ) ) 7 
z{d) — — co7^"^ + 2:07 sec Oq cos{6 - Oq), 



(3.26) 
(3.27) 



where the ray emanates from (0, zq) with initial elevation angle ^Oi and Z07 = ^0 + co7 "'^• 
Arclength is given by s = zq^O, and eigenrays to (r, z) are determined by the angles 



6*0 = 2atan2f2rzo^ - {r^ + {zq-, ~ z)^){r^ + (zq-^ -I- z)^), Zq^ 
6 — atan2(r — zqj tan6'o, z + 007^""^) + Oq, 
with travel-time 



2 z2 



1 



i=: -log . _ 

T \ tan i sin ( 



tan i sin ^ cos 6*0 



^^a±2^cosWo 

Co +7^0 '^ J 



■,(3.28) 
(3.29) 

(3.30) 



assuming no caustics along the ray. Because Cnn = 0, the Hamilton equations (3.3) 



(3.4) are easily solved in closed form, 

g = zo-y sec^ 0o( (1 — a^(co + 7^)^) — |sin6'o|), and p 



1 



Co + 72:0 



(3.31) 
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Likewise, because the Gaussian curvature K = —7^ is constant (the propagation 
medium is a space of constant negative curvature, i.e., it is a hyperboUc space), the Ja- 
cobi equation (3.7) has the closed-form solution f (t) = 7"^ sinh7i, or lit) K,t+ g7^i^ 
for t ^ 7" "'", yielding q = 07^^ sinh7i ss ct{l + 7;7^t^), a much simpler expression than 
Eq. (3.31 ). The expression for the relative delay is given hy p/q w cj~^c^^i^^(l — ^7^^^), 



a new result. 

3.4.3. SSPs with Constant Curvature. The linear SSP example above is 
an example of a medium with constant (negative) Gaussian curvature — a hyperbolic 
space. In general, it is worthwhile to ask which SSPs yield spaces of constant posi- 
tive or negative curvature, because these SSPs serve as models for both convergent 
and divergent ray propagation. Solving the 2d order nonlinear differential equation 
cc" — (c')^ — K for constant K yields, 

cl/W"^ > 0, 



(z) 


— Co cosh(z 


- zo)/W; 


(z) 


— Co sinh(z 


- zq)/W; 


(z) 


— Co cos(z - 


zo)/W; 


(z) 


— Co sin(z — 


zo)/W; 


(^) 


= Co + 7^; 





K 
K 
K 

K = -cl/W"-' 






<0, 
<0, 

<o, 



-7 



<o, 



(3.32) 



zoY is 



of which only the first has constant positive curvature. 

The hyperbolic cosine SSP c{z) = Cocosh(2; — zq)/W « cq + ^coW~'^{z 
the solution for an acoustic duct, which traps rays near depth zq. Regions modeled by 
hyperbolic cosine SSPs are described in the literature with increasing frequency, [27l 
[TllHHIISlIlIlIISllllinilll]; Bergman [5] observes that this profile yields constant 

curvature^ The solution to Jacobi's equation for intrinsic geometric spreading in 

Theorem 



3.4 



equals f (t) ~ K ^ sinK^t, with half-wavelength travel-time between 

caustics determined by ttK~2 — ttW/cq [cf. Eq. (3.12 1]. The range between caustics 
equals co times the travel-time for all small initial elevation angles Oq] therefore, 
by invariance, the range between caustics for all initial elevation angles necessarily 
equals ttW, yielding the following new and useful theorem p71 171 [251 15] : 

Theorem 3.6. Let c{z) = cocosh(z — zo)/W^ be a SSP with hyperbolic cosine 
profile. Then for any initial elevation angle 9q, speed cq, and depth zq, the range 
between caustics and focusing regions is given hy the constant, 



half-wavelength 
range 



TTCo 



R-^^^ = t:W. 



(3.33) 



This theorem is illustrated in Fig. |2.2| over a large range of initial elevation an- 
gles. Porter [3T] uses a value of W~^ = 0.0003 m"^ for his hyperboHc cosine SSP 
and observes without explanation that "the rays refocus perfectly at distances of 
about every 10 km." Indeed, Theorem |3.6| establishes that all rays refocus perfectly 
every 7r/0.0003 = 10.47 km, precisely as shown in Porter's Fig. 8 on p. 2020. [21] 
This is also verified by the closed-form functional form of rays with SSP, available 

from Snell's law, cosh(z — zo)/W^ — sec6'o(l — sin^ 0oCOS^(r/W^)) ^ ; the property of 
constant curvature also yields a simple solutions for the intrinsic geometric spreading 
^{t) = (cq/W) sin(coi/W^) and both the intrinsic phase i^/ ^ = (W/cq) cot{cot/W) and 



extrinsic phase p/g via Eq. (3.22). In practice, it is the quadratic term cosh 2; ^ 1 + ^z"^ 



of the hyperbolic cosine profile that dominates ray behavior in a duct, and ducting 
regions with quadratic SSPs and approximately constant positive curvature are quite 
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common. This establishes the approximate distance between caustics in a duct given 



inEq. (3.12) 



In all other solutions with physical positive speed c(z) > 0, the Gaussian curvature 
is a negative constant, and rays diverge from each other according to the intrinsic geo- 
metric spreading ^{t) ~ A'^3 sinhK^t sa i + ^Kt^, which equals t + ^{cl/W^)-t^ for 
the downward refracting SSP c{z) = cq sinh(z — zq)/W and the trigonometric SSPs 
c{z) = cqCOs{z — zq)/W with divergence zones. The trigonometric SSPs of constant 
negative curvature with convergent ducts (such that c" > 0) are in fact nonphysical 
because the sound speed must be negative at such depths; physical ducting behavior 
with positive sound speed is described by the first hyperbolic cosine solution. The 
case of linear SSPs, whose constant curvature is negative, is treated in the previ- 
ous subsection. Finally, SSPs with vanishing constant curvature K = and intrinsic 
geometric spreading Ht) = t are obviously given by the SSP with constant speed cg. 



3.5. Beam Amplitude and Phase Based Upon Snell's Law. From Fig. |3.2 
it is graphically obvious that the geometric spreading q{t) — ||(9x(t; 0o)/96'o)t||2 sat- 
isfies the trigonometric relationships 



cl- 




dr 



dz 



cos 6, 



(3.34) 



where 9 = tan" ^(i/f) is the elevation angle of the ray at point t. This figure is 
drawn using an implicit assumption Gauss's lemma, namely that the tangent vectors 
X and {dx/d9o)t are perpendicular. It is therefore instructive in establishing the 
relationship between Snell's law and the intrinsic approach to show how Gauss's law 



imply Eqs. (3.34 1. An additional benefit will be an equation for the transverse phase 



of a Gaussian beam expressed in Snell's law form. 

The first-order Taylor series expansion of the expressions r(t; Oq) and z(t; Oq) yield 



Sr^r St^ + {dr/d0o)t SOq, and Sz^ zSt+ {dz/d9o)t SBq. 



(3.35) 



Fixing the variables z and r, i.e., Sz — and 5r ~ respectively, yields the relation- 
ships 



{dr/dea)z = -cot6i • {dz/d9o)r = {dr/ 860)1 - cot 6* • {dz/d0o)t 



(3.36) 



between partial derivatives. Gauss's lemma [Eq. (3.15)] provides the additional rela- 
tionship 



idr/deQ)t + tan6l • idz/d9o)t = 



(3.37) 



that along with Eq. (3.36) yields Eq. (3.34). The conclusion is that the intrinsic 
geometric spreading 



^ = c-\dr/deQ), sin 61 = -c-\dz/d9o)r cost 



(3.38) 



in Snell's law form satisfies Jacobi's equation [Eq. (3.7)], where r{z;0o) is given by 



SneU'slaw [Eq. (2.10)]. 



A benefit of this intrinsic viewpoint of Snell's law is a new expression for the 



transverse phase of a Gaussian beam, obtained by direct application of Eq. (3.22) 
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to Eqs. (2.10) and (3.34|. Differentiating the first equality of Eq. (3.34) w.r.t. s and 



applying the identities dz/ds = sin 9 the Christoffel equations (2.7) yields 

P/l^-\l^] Utt siii^+^^^' (3-39) 



where (dr/dOo)'^ 
viaEq. ( |2.10| : 



{d/dz){dr/ddo)z as usual. These terms are immediately available 
sin 6*0 f^ c{z,) 



dr 
dr 



c(2o) Jz„ (l-a2c2(z,))^/^ 
sin 00 c{z) 



dz,, 



c(2o) (l-a2c2(z)) 



3/2 ■ 



(3.40) 
(3.41) 



Likewise, differentiating the first equality of Eq. (3.38) w.r.t. t yields ^/ 1— c^p/g — c'sin( 
consistent with Eqs. (3.6) and (3.22). 



3.6. Beam Amplitude and Phase Based Upon Ray Angle. The derivation 
of the phase of a Gaussian beam above leads to another formulation of the amplitude 



of a Gaussian beam, which has been used in Gaussian beam applications. Eqs. (3.3) 
and (3.21 ) implies that up to second order in 69q, 



q{s) = / {d9/d0o)s, ds„ 
Jo 



(3.42) 



a formula that has been used by others to compute the extrinsic geometric spreading. 
The transverse phase of a Gaussian beam, up to second-order, is then given by 

cp ( r (d6\ y' foe 



(3.43) 



andEq. (3.18) 



3.7. Intrinsic Transverse Amplitude and Phase of a Gaussian beam. 

The amplitude and phase of a Gaussian beam has been fully developed in Section |3] 
for the computationally practical case of a Gaussian beam defined extrinsically along 
a straight line normal to the ray. In this section, it is explained how Jacobi's equa- 
tion quantifies precisely the same construct intrinsically, if the extrinsic straight line 



in (r, z)-space is replaced with an intrinsic geodesic normal to the ray. Fig. 3.3 il- 



lustrates the difference: because "straight lines" in (r, 2;)-space are actually curved 
intrinsically — they are not time minimizing geodesies — there is a small difference be- 
tween the differential path times. This accounts for the extra term c' sin 6 seen in Eqs. 



(3.6), (3.22), and (3.24). Though computationally impractical, the following new the- 



orem connects the Jacobi equation with the differential path time along rays normal 
to the ray at the center of the Gaussian beam. 

Theorem 3.7. (Intrinsic transverse phase of a Gaussian beam, 2-d case) Let 
Sti be the travel-time difference between the surface of constant travel-time t and an 
intrinsic distance 6fi along the intrinsic transverse plane (i.e., defined by geodesies) 



to the ray at point :ic(t; 9o) , and let ^{t) be the solution to Jacobi's equation [Eq. [3.1)] 
along the ray. Then 



6ti 



14 
2t 



Sn^ + Oidfi^). 



(3.44) 
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Fig. 3.3. Intrinsic view of ray tube in a horizontally stratified medium, in the tangent plane 
at the point (0,2(0)). The azimuthal component, not shown, is perpendicular to the page. All 
distances in the figure represent travel-time, not Euclidean distance. Note that the extri nsic r ay 
normal N = dT/ds — the straight line in extrinsic coordinates used to define Sq in Eg. 1(2.14]) — 
is not straight using intrinsic coordinates and diverges a small amount from the time-minimizing 
geodesic in the normal direction, shown here as the vector Y. The extrinsic normal vector N 
is preferred for numerical computation of the differential path time Stc = \(p / <i) irj^ ; however, the 
intrinsic normal vector Y, from which the differential travel-time is given by Sti = 5 ("^"^f /f ) ^^^ ~ 
\{p/q — c' sin9),Srp , appears in the intrinsic formulation of Gaussian beam. The path time along 
either intrinsic Y or extrinsic N equals Sq = yS6o to first-order. 



The proof of Theorem |3.7| involves another apphcation of Gauss's lemma and a 



Taylor series expansion about x(i;6'o), and is comparable to the proof of Eq. (3.181 
above. Consider the change of elevation angle 9{t) — tan^^(i/r) of the ray: let 69 
be the difference between 9 and the slope of the ray x(i; 6*0 + 69) intersecting the 
transverse plane at distance 6fi = ^69q along a geodesic emanating perpendicularly 
from the ray with initial unit tangent vector Y, i.e., g(Y,Y) = 1 and g{T,Y) = 0. 
This definition of 69 equates to the relation 



(1^, 



'(5p 



±{t;9o + 69)-±{t-9o)), 



(3.45) 



where t7 is the inverse parallel translation along the geodesic between the points 
x(i; 6*0) and x(i; 9o + 69) in direction Y . A straightforward Taylor series analysis 
of the differential equations describing parallelism |13j shows that for any arbitrary 
vector field X{6ii), T^^X{6fi) = X{0)-\-O{6ii'^) — the proof involves so-called "normal 



coordinates" defined so that at the single point where 6/1 = 0, <?ij(0) = 6ij and rays 
in the direction Y have the coordinates (ty^,tF^), implying that Fj (C 
implication of this local analysis about the point x(t; ^q) is that 



0. The 



69^^690 + 0(691) 



(3.46) 



because, up to second order. 



i(i; 9„ + 69) = (d/dt) (x(i; ^o) + tit)Yit)) = i(f ; 9o) + ^Y + ^VtY = i(i; ^o) + ^^ 

(3.47) 
by the parallelism of Y along the ray. Now the change in slope 6ti of the ray is equal to 
the slope of the quadratic function 6ti = {l/2)Qi 611^ + 0{6ijl^) for some coefficient Qi 
that is implied by Gauss's lemma. The slope of this function at (5/i equals 



69 + 0{69^) ^ ta.n 69 = d{6t{)/d(6 ft) ^ QJ^^ Q;^69q. 



(3.48) 



Equating Eqs. (3.46 1 and (3.48 1 , we conclude that Qi = F^/ F, establishing the theorem. 



The quadratic coefficient f /f for the phase is physically consistent with two easily 
imagined examples. For flat space where the Gaussian curvature vanishes, ^{t) = t, 
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Ht) = 1, and the wavefront of constant i is a circle of radius t, in which case the 
gap size equals SU - t sec 60o - t = {t/2){5eaf + 0{h*) = [l / 2){l / 1){5 ^f +0{h'^), 
as predicted by Theorem |3.7[ For the second example, consider acoustic propagation 
in a space of constant positive curvature K = R~^, e.g., S-waves emanating from a 
(geologically unlikely) earthquake on the north pole of the earth. In this case the rays 
are all great circles on the sphere that intersect at the north pole, i.e., longitudinal 
lines. At every point around the equator, equivalently, in the plane transverse to 
every ray at the equator, all rays have equal travel-time to the north pole; therefore. 



Sti = around the equator. Indeed, solving Jacobi's equation [Eq. (3.7)] yields 
lit) — sin{t/R) for which f' = cos(|7r) = at the equator, an easily obtained intuitive 
result that is again consistent with Theorem |3.7[ 



4. Conclusions. Paraxial ray theory is an essential and useful component of 
modern acoustic modeling. There are deep and important connections between the 
classically derived paraxial ray equations for the amplitude and phase along a Gaus- 
sian beam, and the second-order variation — called the Jacobi equation — along geodesies 
encountered in differential geometry. This paper demonstrates how known results in 
paraxial ray theory correspond to their counterparts in differential geometry, and 
shows how both new equations and new insights into the properties of acoustic rays 
are obtained from an intrinsic, differential geometric point of view. It is shown how 
the intrinsic Gaussian curvature affect the spreading of Gaussian beams, and how the 
specific form of this intrinsic curvature [K = cc" — (c')^ for a horizontally stratified 
medium with SSP c(z)] allows one to easily compute the distance between caustics 
within a duct, as well as the geometric spreading for either a duct or a region with 
linear SSP. These results allow the introduction of SSPs yielding constant Gaussian 
curvature, which serve as model spaces for both convergent acoustic ducts (positive 
curvature), divergence zones (negative curvature), and non-refractive spreading (zero 
curvature). It is proved for the model of a hyperbolic cosine SSP, the range between 
caustics is constant for all rays. Intrinsic versions of the amplitude and Gaussian 
beams are introduced. The intrinsic geometric spreading is shown to be equivalent 
to its extrinsic counterpart after scaling by the position-dependent sound speed. The 
intrinsic phase of a Gaussian beam, which is the phase along geodesies, not straight 
lines, emanating in a normal direction from the ray, is shown to be equivalent up to a 
small additive term to the phase of a paraxial ray as it is typically defined. Because 
the differential geometric approach is quite general, all results may be generalized to 
three dimensions, and in all cases, the connection is made with known results de- 
rived from cither Snell's law, Hamilton's equations, or the Frenet-Serret formulae. 
The results may also be applied to other special cases of applied interest, such as a 
spherically stratified sound speed. 
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